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Abstract: We consider the possibility that the muon g — 2 anomaly, S {g — 2)^, finds its 
origins in low energy supersymmetry (SUSY). In the general MSSM the parameter space 
consistent with S {g — 2)^ and correct dark matter relic density of the lightest neutralino 
easily evades the present direct LHC limits on sparticle masses and also lies to a large 
extent beyond future LHC sensitivity. The situation is quite different in GUT-defined 
scenarios where input SUSY parameters are no longer independent. We analyze to what 
extent the LHC can probe a broad class of GUT-inspired SUSY models with gaugino non¬ 
universality that are currently in agreement with the bounds from 5 {g — 2)^, as well as 
with the relic density and the Higgs mass measurement. To this end we perform a detailed 
numerical simulation of several searches for electroweakino and slepton production at the 
LHC and derive projections for the LHC 14 TeV run. We show that, within GUT-scale 
SUSY there is still plenty of room for the explanation of the muon anomaly, although 
the current LHC data already imply strong limits on the parameter space consistent with 
S (g — 2)^. On the other hand, we demonstrate that the parameter space will be basically 
fully explored within the sensitivity of the 14 TeV run with 300 fb“^. This opens up the 
interesting possibility that, if the {g — 2)^ anomaly is real then some positive signals must 
be detected at the LHC, or else these models will be essentially ruled out. Finally, we 
identify the few surviving spectra that will provide a challenge for detection at the LHC 
14 TeV run and we characterize their properties. 


^On leave of absence from the University of Sheffield, U.K. 
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1 Introduction 

The measurement of the anomalous magnetic moment of the muon, [g — 2)^, by the 
Brookhaven experiment [1] a decade ago shows a ~ 3cr discrepancy with the Standard 

Model (SM). The measured value of this discrepancy, S {g — 2)^ = — o™, is 

(5(5-2)^ = (28.7 ±8.0) X 10"^° or (1.1) 

(5(5-2)^ = (26.1 ±8.0) X 10"^°, (1.2) 

according to whether the lowest order hadronic contributions from Ref. [2] or Ref. [3] are 
used to compute the SM value. The discrepancy will soon be either confirmed or overruled 
by the New Muon g-2 experiment at Fermilab [4, 5] , which is bound to rekindle the interest 
of the particle physics community in the subject. 

A ~ 3(T difference with the SM can easily be accommodated in the Minimal Supersym¬ 
metric Standard Model (MSSM) [6-32], which is rich in particles with the right couplings 
to provide significant loop contributions to the ggL'y vertex. The 1-loop contributions in 
the MSSM [7, 33] are roughly split into those arising from a chargino/sneutrino loop and 
those given by smuon/neutralino loops, so that at the lowest order the measurement of 
6 {g — 2)^ can be parametrized by 

p,Mi,M2,m/i^,m/i^,m£;^,tan/3, (1.3) 

where g is the higgsino mass parameter, tan /3 is the ratio of the Higgs vacuum expectation 
values (vev’s), Mi and M 2 are the bino and wino soft supersymmetry-breaking masses, 
and are smuon soft masses, and is the soft mass of the muon sneutrino. 
The LHC has started to test the electroweak (EW) sector of the MSSM in a class 
of searches involving different multiplicities of leptons, no jets, and a significant amount 
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of missing energy. The 8 TeV run has provided strong lower bounds on chargino and 
slepton masses, in particular when interpreted in the framework of simplified model spectra 
(SMS) [34], However, several studies have shown [19, 22, 23, 28, 29] that if the experimental 
limits provided by the CMS and ATLAS collaborations are reinterpreted and applied to 
more general MSSM scenarios the 8 TeV LHC results can only constrain a small part of the 
available parameter space, so that ample room still remains to attribute a supersymmetric 
(SUSY) origin to 5 (g — 2)^. 

On the other hand, in scenarios where unification of the scalar and gaugino masses is 
imposed as a remnant of some new physics at the scale of Grand Unification (GUT scale), 
like in the well-known Constrained MSSM (CMSSM) or the Non-Universal Higgs Mass 
(NUHM) model, it has become virtually impossible to find regions of the parameter space 
where the measurement of <5 (<7 — 2)^ can be reproduced (see, e.g., [35-38]). The reason is 
well known, and is due to the combined effect of direct lower limits on colored sparticles at 
the LHC and the discovery of a Higgs boson with nih — 125 GeV [39-41], which together 
have pushed the favored parameter space for unified scalar and gaugino masses to the 
multi-TeV regime, thus indirectly forbidding the possibility of light sleptons, binos, and 
winos. 

Interestingly, this is not necessarily the case for GUT-scale models in which the as¬ 
sumption of gaugino unification is relaxed. It has been pointed out, e.g., in Refs. [22, 42] 
(but see also [20, 27, 31, 43]), that in models of gravity mediation all one really needs is 
GUT-defined boundary conditions such that the high-scale value of the gluino soft mass, 
M 3 , is much larger than the values of Mi and M 2 . Through the renormalization group 
equations (RGBs) large initial conditions for M 3 can then drive the masses of sparticles 
charged under SU{3) to large values at the EW scale, in agreement with the LHC data 
and the measurement of the Higgs boson, while the sparticles charged under only the EW 
gauge groups remain quite light, in agreement with the measurement of 6 {g — 2 )^. 

The goal of this paper is to examine to what extent the oncoming run of the LHC 
can probe the parameter space of these gravity-mediated, GUT-scale SUSY models that 
satisfy the present constraints for {g — 2)^. Besides being motivated by considerations of 
symmetry and providing a realistic framework for SUSY breaking, these scenarios are more 
constrained than generic phenomenological parametrizations of the MSSM by the measure¬ 
ment of the relic abundance of dark matter (DM) at Planck or WMAP [44, 45], VLh? ~ 0.12. 
Once this bound is combined with the measurement of 5 {g — 2 )^, the allowed parameter 
space becomes significantly reduced. Eor example, we shall see that often the lightest 
SUSY particle (LSP) is a fairly light bino-like neutralino y, which needs an equally light 
slepton or, alternatively, substantial mixing with higgsinos to yield < 0.12. Thus, 

the parameter space that gives and 6 {g — 2 )^ is in this case particularly sensitive to 

direct LHC searches involving the production and decay of sleptons and electroweakinos. 

In this paper we reinterpret existing 3-lepton [46, 47] and 2-lepton [48] LHC searches 
for direct production of charginos, neutralinos, and sleptons, and apply them to the allowed 
parameter space of several GUT-scale SUSY models characterized by non-universal bound¬ 
ary conditions to fit the {g — 2)^ measurement. We progressively increase the complexity 
of the analyzed models by disunifying several parameters defined at the GUT scale. We do 
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this by following patterns typical of GUT symmetries like ^©(lO), SU{5) or Pati-Salam. 
From the phenomenological point of view, this is equivalent to introducing an increasing 
number of mechanisms that yield the correct relic abundance of the LSP. 

We will show that the present constraints from the LHC on the EW sector of this 
large class of models are already quite stringent and, more importantly, that the 14TeV 
run offers the opportunity to probe the remaining parameter space virtually in its entirety. 
To this end, we calculate the projected sensitivity of 2- and 3-lepton searches for the LHC 
14TeV run with 300 fb~^ and apply the simulations to the defined GUT-scale models. 
Incidentally, our projections can be compared for SMS scenarios to the ones produced by 
other groups [49]. 

The paper is organized as follows. In Sec. 2 we review the interplay of the bounds from 
the measurement oi 6 {g — 2)^ and on the parameter space of the MSSM, and we will 
comment on the LHC reach in the phenomenological MSSM. In Sec. 3 we introduce the 
GUT-scale models that will be scanned over and expose the parameter space consistent 
with several phenomenological constraints. In Sec. 4 we present our methodology for 
reinterpretation of LHC searches and our projections for the next run of the LHC in SMS 
scenarios. We then show the main results, i.e., we apply the derived LHC bounds and 
projections to the GUT scenarios defined in Sec. 3. We finally present our summary and 
conclusions in Sec. 5. 


2 The relic density and (g — 2)^ in the MSSM 


We review in this section the relations that can be derived on the parameter space of the 
MSSM when the measurement of 6 (g — 2)^ is considered in combination with the relic 
density. A study on the subject that includes the limits from the LHC 8 TeV run is done, 
e.g., in Ref. [28]. Our goal here is to show that these relations can be very useful in deriving 
bounds, but unless some assumption about the mechanism of SUSY breaking is made, there 
remain large fractions of the parameter space outside of the LHC reach. 

The MSSM contributions to 5 {g — 2)^ have been calculated at 1 loop in [7, 33, 50- 
53] using the mass insertion method.^ Higher order contributions have been computed 
in [55-63]. At 1 loop there are five main contributions that can be split into two classes of 
diagrams: chargino/sneutrino and neutralino/smuon contributions. In terms of the MSSM 
parameters these five contributions are given by [33]: 


tan/3 / M| \ 

(47j-)2 gM2 J ’ 

(1) ^ _l^^m^tan^ _ / M| \ 

2 (dvr)2 fiM 2 J ’ 

a( 2 ) ^ 1 g'^ m^tan/3 _ / Mf \ 

2 (47r)2 g,Mi ' 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


code to calculate the 1-loop contributions to 5 {g — 2) in generic new-physics models has recently 


become available [54]. 
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Table 1: Regions in the MSSM giving 5 {g — 2)^ (1-loop) at 2a and < 0.12. The symbol > means 
here “greater but not orders of magnitude greater than...,” see Figs. (l)-(3). 
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where 5 and g' are the gauge couplings of the SU{2) and (7(1) SM groups, respectively, 
and the j and are loop functions that read 


T, 
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( 2 . 6 ) 

(2.7) 


where we have used the reduced forms of Ref. [19]. Note that the numerical coefficient 
in front of Eqs. (2.3) and (2.4) depends on g'^ so that these contributions are in general 
suppressed with respect to Eq. (2.1). The neutralino/smuon contribution of Eq. (2.5), 
however, depends directly on //. When g it can become the dominant 

one. 


The bounds on the parameters (1.3) that arise from the measurement oi 5 {g — 2)^ can 
be combined with the bounds that come from imposing < 0.12 on the same regions 

of the parameter space. We review here the regions consistent with Eq. (1.1) at least at 
the 2 (T level and show the correspondingly allowed parameter space. 

We adopt in this section simplifying assumptions typical of many phenomenological 
parametrizations of the MSSM [64]. The soft SUSY-breaking parameters are defined at the 
EW scale and we assume that first and second generation slepton soft masses are degenerate 
(mg^ = mjx^ = and similar identities apply to right-handed sleptons). 

The bounds obtained using the 1-loop calculation, Eqs. (2.1)-(2.5), are not far off from 
the ones given by more precise higher-order calculations, so that for the semi-quantitative 
discussion of this section we will limit ourselves to the former approximation. In the 
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numerical analysis that we present in the following sections we will consider higher-order 
contributions with the help of the latest numerical codes. 

The bounds arising from the combination of {g — 2)^ and are summarized in 

Table 1 and discussed below. 


.Z//i-resonance. As is well known, when the lightest neutralino (hereafter simply “the 
neutralino” or x) is in the mass range 30 GeV < < 62 GeV the cross section for pair 

annihilation can be enhanced by the resonance with the Z boson or the Higgs: XX 
Z/h —7- SMSM [65, 66]. The neutralino is in this case predominantly bino-like. 

To undergo Z-resonance annihilation, when ~ M 2 / 2 , the neutralino must have a 
non-negligible higgsino component to maximize the coupling to the Z boson, so that the 
cross section is given by, up to proportionality constants and phase-space integration, [65, 
67] 


ji 


(TV 


mf 


1 - 


M? 


4- 
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FzMz^ 


2 ’ 


( 2 . 8 ) 




where Tz — 2.5 GeV is the Z width. Obviously, when /i becomes large the cross section 
decreases and the relic density exceeds the measured value. This effectively sets an upper 
bound, g < 500 GeV corresponding to ~ 0.12. When M 2 > /i one derives an upper 

bound on the mass of the (mostly higgsino-like) lightest chargino, m + < 500 GeV. 

Xl 


The Higgs resonance [66] is qualitatively similar to the Z-resonance with a few differ¬ 
ences. The width of the Higgs boson is much narrower than the Z’s: T/j ~ 4MeV [68], so 
that the cross section is more sensitive to the neutralino being on or off the resonance. But, 
more importantly, the cross section scales as 1/(1 —/r/Mi)^ [67] rather than 1/{1 — g? 
as in Eq. (2.8), so that g is less constrained than in the Z-resonance region, g < 1000 GeV 
in the /i-resonance region. On the other hand, Z- and /i-resonance regions both depend 
minimally on tan /3 or the slepton masses, so that these parameters are not bounded by 
the relic density constraint. The same is true for the wino soft mass, M 2 , that can assume 
arbitrarily large values. 

Because of the upper bound on g, in the Z and /i-resonance regions the dominant 
contribution to 6 {g — 2)^ is given by the chargino-sneutrino loop, Eq. (2.1). The parameter 
space allowed at 2a by Eq. (1.1) strongly depends on the value of tan/3, which cannot be 
constrained by 

In Fig. 1(a) we show in the {m ±, rrij ) plane the part of the Z//i-resonance region 
consistent at 2a with Eq. (1.1). The colored bands show the {g — 2)^ bounds for two values 
of tan 13, tan /3 = 60 in orange and tan /3 = 10 in yellow. The right-handed slepton mass 
is set to a large value, = 5 TeV, as the dominant contribution to 6 {g — 2)^ in this 
case, Eq. (2.1), does not depend on fni^. The plot shows that soft mass is bound to be 
lighter than ~ 1.5 TeV for tan /3 ~ 60, but the upper bound becomes more restrictive as one 
considers smaller tan j3 values. We assume here M 2 = 2g. For larger values of M 2 the plot 
gets slowly squashed down, and at M 2 — 5 TeV the limits on become approximately 3 
times stronger. 

It has been shown [24, 28, 69] that the parameter space corresponding to the Z/h- 
resonance region can be probed at the LHG by 3-lepton searches for EW-ino production 
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Figure 1: (a) The parameter space belonging to the Z//i-resonance region of the MSSM that is allowed 
at 2 (t by (g — 2)^ in the ijn ^±, rri;-^) plane. The orange band shows the case with tan /J = 60 and the yellow 
band the one with tan j3 = 10. The black dot-dashed horizontal line shows the approximate 95% C.L. lower 
bound for from 2 lepton searches [47, 48] at the LHC. The red dashed line shows the approximate 
lower bound in from 3 lepton searches [46, 47], which differs if or viceversa. (b) The 

allowed 2a parameter space in the ) plane for the parameter space belonging to the region of mixed 

bino/higgsino neutralinos. The solid dark blue line shows the approximate 95% C.L. exclusion bound from 
2-lepton searches. The color code is the same as in (a). 
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and 2-lepton searches for direct slepton production. The bounds from 3-lepton searches are 
much stronger in the presence of an intermediate slepton between the mass of the chargino 
and neutralino [46, 47]. They are approximately shown as red dashed lines in Fig. 1(a). The 
approximate bound from 2-lepton searches [47, 48] for slepton pair production is shown as 
a dot-dashed black line. One can see that for large tan (3 a large fraction of the parameter 
space is presently not excluded. However, we will show in Sec. 4 that the outlook for the 
14 TeV run improves considerably. 

Neutralino of mixed bino/higgsino composition. As one considers larger masses, for 
a bino-like neutralino it becomes necessary to increase the mixing with higgsino states to 
enhance the annihilation cross section. These “well tempered” neutralinos [70] efficiently 
annihilate to gauge bosons through t-channel chargino exchange. 

In Fig. 1(b) we show in the {'^ 1 ^, my.) plane the region of the parameter space consistent 
at 2a with Eq. (1.1) in the case of a mixed bino/higgsino neutralino with rriy. ~ Mi < //. 
The colored bands show two different tan/3 cases. Again the right-handed slepton mass 
has been set at 5 TeV, and the wino soft mass is set here to M 2 = Mi -\- 500 GeV. Raising 
M 2 moves the allowed bounds down and left, to smaller values of rriy. and by reducing 
the contribution of Eq. (2.1) to <5 (<7 — 2)^. Thus, M 2 cannot be heavier than a few TeV. 
The {g — 2)^ constraint requires approximately, < 1.2 — 1.5 TeV and < 700 GeV. 
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This region is notoriously difficult to probe in 3-lepton searches at the LHC, because 
the masses of the lightest chargino and neutralino are almost degenerate. On the other 
hand, the bounds from 2-lepton searches, approximately indicated with a solid dark blue 
line, are much too weak at the moment to probe the parameter space. 

Neutralino/slepton coannihilation. For a predominantly bino-like neutralino the cor¬ 
rect relic density can be obtained if y coannihilates with an almost degenerate slep- 
ton [71, 72]. This mechanism is particularly important for what follows, as it is one of 
the few realized in models with GUT-scale boundary conditions. 

In most cases the neutralino mass must be within ~ 20 GeV of the mass of the coan¬ 
nihilating slepton, and it also becomes very hard to compensate for the increasing mass 
when rriy' reaches an approximate upper bound of ~ 500 GeV. 

In Fig. 2(a) we show the 2 ct allowed parameter space in the m^) plane for 

the stau-coannihilation region, in which the neutralino coannihilates in the early Universe 
with the lightest stau. The dominant contribution to S {g — 2)^ is given in this case by 
the neutralino/smuon loop of Eq. (2.5), which increases linearly with g. Thus, values of 
g much larger than in the previous cases are allowed and they actually help to satisfy the 
{g — 2)^ constraint. The bounds for two very different values of g are shown in the plot, 
g = I TeV (yellow band) and g = lOTeV (orange band), while tan/3 is kept large. 

In Fig. 2(b), the value of g is instead fixed at an intermediate value, g = 5 TeV, and we 
show the bounds for two different values of tan (3: tan /3 = 10 (yellow band) and tan /3 = 50 
(orange band). 

We show in both panels the case with M 2 = 5 TeV and = rrij -^, to maximize the 
contribution of Eq. (2.5). Note that since Eq. (2.5) does not depend on M 2 , the wino soft 
mass can be actually decoupled and without additional assumptions on the mechanism 
of SUSY breaking 3-lepton searches are not in principle sensitive to this region of the 
parameter space. On the other hand, one can see in Eigs. 2(a) and 2(b) that the {g — 2)^ 
constraint bounds the value of the smuon (and selectron) masses, but the sensitivity in 
2-lepton searches is still very limited to bite significantly into the parameter space. Note, 
finally, that the {g — 2)^ bounds on the left- and right- handed slepton masses presented 
in Figs. 2(a) and 2(b) become weaker in cases where g and M 2 are both beneath ~ I TeV, 
as Eq. (2.1) becomes then dominant. 

Rather than with the lightest stau, the neutralino can coannihilate with a light se¬ 
lectron, smuon, or sneutrino. The {g — 2)^ bounds in the 'nrix) pis'll® for the case 

of coannihilation with a left-handed slepton of the first or second generation is shown in 
Eigs. 2(c) and 2(d). One can obtain similar plots, which well approximate the case of coan¬ 
nihilation with the right-handed selectron or smuon, by replacing —)• (Eq. (2.5) 

is symmetric under •(-)> ?7i[^)- Again we show in Fig. 2(c) the case of large tan/3 for 
different values of g, and in Fig. 2(d) the case of fixed g for different values of tan/3. For 
coannihilation to occur the left-handed slepton mass is kept relatively low, ~ m^, so 
that the (^f — 2)^ upper bounds on the right-handed mass are actually weaker than in the 
stau-coannihilation region: ^ 2 — 2.5 TeV in this case. 

The dark blue solid lines in Fig. 2 show the approximate bound on the left-handed 
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Figure 2: (a) The parameter space allowed at 2a by {g — 2)^ in the m^) plane for the stau- 

coannihilation region of the MSSM. We assume a large value of tan /3. The orange band shows the case with 
/i = 10 TeV and the yellow band the one with /r = 1 TeV. The solid dark blue line shows the approximate 
95% C.L. exclusion bound from 2-lepton searches at the LHC. (b) Same as (a), for a fixed value /r = 5 TeV. 
The orange band shows the case with tan /3 = 50 and the yellow band the one with tan /3 = 10. (c) The 
parameter space allowed at 2a by {g — 2)^ in the mx) plane for the left slepton-coannihilation region 

of the MSSM with fixed tan/3 = 50. The color code is the same as in (a), (d) The parameter space allowed 
at 2a by {g — 2)^ in the plane for the left slepton-coannihilation region of the MSSM with fixed 

= 5 TeV. The color code is the same as in (b) . 


selectron mass from the LHC 2-lepton searches for pair production. In the cases 

shown in Figs. 2(c) and 2(d) the reader should take note of a couple of caveats: the first 
is that the bound does not properly apply when a left-handed selectron is degenerate with 
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the neutralino because the spectrum is compressed; the second is that the limit for right- 
handed slepton pair production is actually weaker than the one shown here by ~ 50%, 
because the cross section for eRCR production is suppressed with respect to left-handed 
production. The proper and complete treatment of the LHC limits for the cases with 
coannihilation will be presented in Sec. 4. 

Nearly pure higgsinos and winos. We finally discuss the parameter space corresponding 
to an apparent underabundance of DM in the Universe, Q.y.h? < 0.12. As is well known, this 
situation is typical when the LSP is an almost pure higgsino with <C 1 TeV [73, 74] or 
an almost pure wino with <C 2.8 TeV [75]. These solutions are generally thought to be 
proper of scenarios with two-component DM, or of cases where the neutralino represents 
the entirety of the DM and the correct abundance is fixed by invoking some additional 
mechanism, e.g., freeze-in [76]. 

We show in Fig. 3(a) the 2(t allowed parameter space in the plane for a 

higgsino LSP. The plot is generated under the assumptions that M 2 is greater than /r, but it 
cannot be orders of magnitude greater, otherwise the dominant contribution to 5 {g — 2)^, 
Eq. (2.1), becomes drastically suppressed, as was the case for the mixed bino-higgsino 
scenario. On the other hand, Eq. (2.1) is insensitive to the values of and Mi, which 
can be decoupled. 

In Fig. 3(b) we show the allowed parameter space for a wino-like neutralino with 
relatively large fi. Besides Eq. (2.1), substantial contribution to 5 {g — 2)^ comes in this 
case from Eq. (2.5), so that 5 {g — 2)^ can be enhanced for large jjL when the slepton and 
bino masses are not much above M 2 ■ For smaller values of fi the behaviour becomes similar 
to that of the pure higgsino and mixed bino-higgsino cases. 

As was the case for the bino/higgsino admixtures described above, almost pure higgsi¬ 
nos and winos are extremely difficult to test at the LHC in 3-lepton final state searches, 
because of the strong degeneracy between and y. On the other hand, it is also clear 
from Fig. 3 that 2-lepton searches can begin to test the parameter space in these cases 
but, especially for larger values of tan/3, it will be very hard to reach enough sensitivity to 
probe the full allowed parameter space. 

To summarize, this discussion has proven that, when one requires consistency with the 
{g — 2)^ constraint and the relic density, important regions of the MSSM parameter space 
are within the reach of the LHC. However, in the MSSM there is too much freedom and one 
cannot entirely test the hypothesis of a SUSY origin to the {g — 2)^ anomaly at the LHC, 
especially in cases where degeneracy among the main parameters becomes important. In 
what follows, we will focus on theoretically well-motivated scenarios where the presence of 
additional symmetries at the GUT scale removes much of the freedom of the parameter 
space. We limit ourselves to the most common case of a bino-like neutralino, which can 
saturate the relic density and is testable at the LHC, but we will comment on the cases 
that can give rise to neutralinos of a different composition. 
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Figure 3: The parameter space allowed by {g — 2)^ at 2a in the m^) plane for some cases with 

< 0.12. (a) Case of a higgsino-like neutralino with ~ p < 1 TeV. (b) Case of a wino-like neutralino 
with nix ~ -^2 ^ 3TeV. In both plots the orange band shows the case with tan/3 = 60 and the yellow 
band the one with tan/3 = 10. The solid dark blue line approximately shows the 95% C.L. exclusion bound 
from 2-lepton searches at the LHC. 
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3 GUT-defined models and experimental constraints 

Many models of gravity mediation assume for simplicity universal GUT-scale conditions 
for the gaugino and scalar soft masses. As was mentioned in Sec. 1, however, under those 
assumptions it has become no longer possible to find parameter space consistent with the 
{g — 2)fj, constraint after the Higgs discovery and null searches for squarks and gluinos at 
LHC run I. As a matter of fact, in scans of the CMSSM and the NUHM the {g — 2)^ 
constraint is in some cases neglected [36, 77] when looking for the regions of the parameter 
space favored by the Higgs measurement or the LHC. The issue is somewhat controversial 
and judgement is in general postponed to after an eventual confirmation by more precise 
upcoming experiments. 

We follow here a different approach and consider GUT-defined models that actually 
do satisfy the present constraints for {g — 2)^. To this end, we relax the assumption of 
gaugino universality. We will show in the remainder of this paper that the LHC bounds on 
the EW sector of these GUT-defined SUSY models are quite strong and that the 14 TeV 
run will be able to probe the parameter space of these models virtually in its entirety. 

The models we analyze are summarized in Table 2 and discussed below. The scans 
are performed with the package BayesFITS [23, 36, 77, 78] which interfaces several pub¬ 
licly available tools to direct the scanning procedure and calculate physical observables. 

The sampling is performed by MultiNest [79] with 4,000 or 6,000 live points. We use 
SoftSusy v.3.5.2 [80] to calculate the mass spectrum and SUSY-HIT [81] for the decay 
branching ratios. Higher-order corrections to the Higgs mass are calculated with FeynHiggs v.2.10.2 
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Model 1 

CMSSM-like Ms floating 


Parameter 

Description 

Range 

mo 

Universal scalar mass 

100, 4000 

"ll/2 

Bino/wino soft mass 

100, 4000 

M3 

Gluino soft mass 

700, 10000 

^0 

Universal trilinear coupling 

- 8000, 8000 

tan /? 

Ratio of the Higgs vevs 

2, 62 

sgn/r 

Sign of the Higgs/higgsino mass parameter 

+ 1 

Model 2 

Non-universal gaugino masses 


Ml 

Bino soft mass 

- 4000, 4000 

Ms 

Wino soft mass 

- 4000, 4000 

mo, Ms, Aq, tan/3, sgn/i 

Same as Model 1 

Same as Model 1 

Model 3 

S'O(10)-like sfermions 


mi6 

Universal scalar mass 16 repr. 

100, 4000 

mfo 

Universal scalar mass 10 repr. 

- 100002, 100002 

2, Ml 

D-term extra U{1) 

0, mfg — (100GeV)2 

mi/2 

Bino/wino soft mass 

100, 2000 

Ms 

Gluino soft mass 

800, 5000 

Aq, tan/3, sgn g 

Same as Model 1 

Same as Model 1 

Model 4 

S'Ll(5)-like sfermions 


mio 

Universal scalar mass 10 repr. 

100, 4000 

ms 

Universal scalar mass 5 repr. 

100, 2000 


Down Higgs doublet soft mass 

- 100002, 100002 

9 

mrj 

Up Higgs doublet soft mass 

- 100002, 100002 

mi/ 2 , Ms, Aq, tan/3 sgn// 

Same as Model 3 

Same as Model 3 

Model 4-zoom 

SU{5) g,,mA parameterization 



EW-scale higgsino mass parameter 

10, 2000 

ruA 

Pseudoscalar pole mass 

100, 4000 

Ms 

Gluino soft mass 

500, 2000 

mio, ms, mi/2, Ao, tan^ 

Same as Model 4 

Same as Model 4 


Table 2: Parameters of the models analyzed in this work. All soft SUSY-breaking masses 
are defined at the GUT scale. Dimensionful quantities are given in GeV and GeV^. 


[82-86]. FeynHiggs is interfaced with HiggsSignals vl.3.1 [87] and HiggsBounds v4.2.0 [88- 
90] to evaluate the constraints on the Higgs sector. SuperlSO v.3.4 [91] is used to calculate 
6{g — 2)^ and flavor observables BR(B—)-Xs 7 ), BR (Bg —)•and BR (B^ —)■ ri^). 
M\y, sin^^effj are calculated with FeynHiggs. Dark matter observables, and 

the spin-independent DM-proton cross section, are computed with micrOMEGAs v.4.1.5 [92] . 
SuperlSO v.3.4 performs the calculation of 6 {g — 2)^ including the leading-log QED 
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Constraint 

Mean 

Exp. Error 

Th. Error 

Ref. 

Higgs sector 

See text. 

See text. 

See text. 

[87-90] 

LUX 

See [77, 93]. 

See [77, 93]. 

See [77, 93]. 

[94] 


0.1199 

0.0027 

10 % 

[44] 

6{g-2)^ X IQio 

28.7 

8.0 

3.0 

[1, 95] 

sin^ 6»eff 

0.23155 

0.00015 

0.00015 

[96] 

BR (B ^ Xs 7 ) X 10^ 

3.43 

0.22 

0.21 

[97] 

BR (Bu — )■ Tu) X 10^ 

0.72 

0.27 

0.38 

[98] 

AMb^ 

17.719 ps-i 

0.043 ps-i 

2.400 ps"^ 

[96] 

M\y 

80.385 GeV 

0.015 GeV 

0.015 GeV 

[96] 

BR(Bs — )■ X 10 ® 

2.9 

0.7 

10 % 

[99, 100] 

r(^ ^ xx) 

< 1.7 MeV 

0.3 

- 

[ 101 ] 


Table 3: The experimental constraints applied in this study. 


corrections from 2-loop evaluations [55], photonic Barr-Zee diagrams with physical Higgs [56- 
58] and bosonic EW 2-loop contributions [60]. 

The scans are subject to a set of constraints, applied through a global likelihood func¬ 
tion C. The list of constraints, central values, theoretical and experimental uncertainties 
are presented in Table 3. We assume Gaussian distributions for the constraints, with the 
exception of those on the Higgs sector, which are imposed through HiggsSignals and 
HiggsBounds, and the constraints on from LUX [94]. The LUX constraint, which 
slightly improved on the limit from XENONIOO [102], is included in the likelihood function 
following the procedure detailed in [23, 93, 103]. Additionally, we impose 95% C.L. lower 
bounds from direct searches at LEP [101], smeared with 5% theoretical errors. The limits 
are given in Eq. (2) of Ref. [23], with the exception of the limit on the neutralino mass that 
has been replaced here by the LEP limit on the invisible Z width, V{Z —)■ xx) [101]. 

To define the 2a allowed regions, we adopt for S {g — 2)^ the central value of Eq. (1.1). 
We estimate the theoretical uncertainty very conservatively, = 3.0 x 10“^*^, to bundle 
together the uncertainties that arise from neglecting hadronic 2-loop corrections [62, 63] in 
the SUSY calculation and the SM uncertainties that give rise to different estimates, like 
Eq. (1.2). 

The first model we consider. Model 1 hereafter, is a simple modification to the 
CMSSM first introduced in Ref. [22].^ In addition to the usual parameters, the gluino 
soft mass M 3 is allowed to float at the GUT scale, as shown in Table 2. This is the 
minimal implementation of non-universality that allows one to simultaneously respect the 
bounds on the color sector from the LHC and the Higgs mass, and those on the EW sector 
from [g — 2 )^. As was shown in Ref. [ 22 ], this simple condition can be easily obtained 
within several GUT symmetries. 

^The authors of Ref. [22] call these scenarios ySUGRA, as the radiative breaking of EW symmetry is 
driven by the gluino, g. 
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The distributions of the input parameters of Model 1 after applying the constraints of 
Table 3 are not particularly illuminating for the purposes of this paper and we refrain from 
showing them here. As expected, at the GUT scale M 3 assumes large values, M 3 ~ 1 — 
5 TeV, constrained by the fact that it must drive the physical stop masses to the multi-TeV 
regime to comfortably fit the Higgs mass and rates. One finds the following approximate 
relations between the physical masses and the GUT-scale value of M 3 : ~ 1.5 M 3 , 

rrig ~ 2 M 3 . Gonversely, the common scalar mass remains small, mo < 450 GeV, due to the 
fact that the physical smuon mass must be small enough to be consistent with S {g — 2 )^. 

As is often the case in global SUSY analyses, the relic density provides the constraint 
with the strongest impact on the EW sector. In Fig. 4(a) we plot with magenta triangles 
the distribution of the physical left-handed selectron mass, nig^, versus the neutralino mass, 
m^, for the points of Model 1 satisfying the constraints of Table 3 at the 2a level. The 
right-handed selectron mass distribution, mg^ is shown with blue circles, and the lightest 
stau mass distribution, m^^, with cyan diamonds. One should remember that in gravity- 
mediated models m^^ ~ given the very small value of the Yukawa couplings, and the 
same is true for the right-handed sleptons. 

The neutralino is strongly bino-dominated and the correct relic abundance is obtained, 
for 100 GeV < < 350 GeV, through coannihilation with the lightest stau. For a lighter 

m^ it is obtained through bulk-like annihilation [65] to taus via t-channel exchange of 
the moderately light stau. One can see that coannihilation rapidly loses efficiency as m^ 
increases, so that no solutions are found for > 350 GeV. 

As was explained in Sec. 2 the dominant contribution to 5 {g — 2 )^ is given by the 
neutralino/smuon loop of Eq. (2.5) so that the parameter /i adopts fairly large values, 
// > 2 — 5 TeV. Unlike in general low-scale MSSM scenarios like the one described in Sec. 2, 
in Model 1 tan/3 can only assume moderate values, tan/3 ~ 5 — 25, as for larger tan/3 the 
stau masses run to unphysical low-scale values giving ruf^ < m^. 

The GUT-scale universality condition Mi = M 2 = mi /2 is relaxed in Model 2, whose 
parameters are shown in Table 2. We show in Fig. 4(b) the distribution of the physical left- 
handed selectron mass, right-handed selectron mass, and stau mass versus the neutralino 
mass for the 2a allowed points. As can be inferred by the stau mass distribution, the 
mechanism of interest for the relic density is again stau-coannihilation, like in Model 1. 
The main difference with the previous case is that, in addition to a broader range of m ± 
values, there is a broader distribution for the slepton masses of the first two generations. 

This can be understood by looking at Fig. 4(c), where we show with magenta triangles 
the distribution of the left-handed selectron mass versus the GUT-scale ratio M 2 /M 1 , 
and with blue circles the corresponding distribution for the right-handed selectron mass. 
Note that we can only find solutions in the 2cr region of the constraints of Table 3 when 
M 2 /Ml > 0. Model 1 is the subset of Model 2 represented by the points at M 2 /Mi = 1. One 
can see that the splitting between the left- and right-handed selecton masses increases with 
increasing M 2 /Mi, as larger M 2 values can drive the left-handed mass to larger values at 
the low scale through the RGEs. Moreover, for ratios larger than M 2 /M 1 ~ 3.5 the fi soft 
mass becomes too large after running to the low-scale to allow for efficient coannihilation 
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Figure 4: (a) The distribution of the physical left-handed selection mass (magenta triangles), right- 
handed selection mass (blue circles), and lightest stau mass (cyan diamonds) versus the neutralino mass for 
the points of Model 1 that satisfy the constraints of Table 3 at the 2tj level, (b) Same as (a) for Model 2. (c) 
The distribution of the physical left-handed selection mass (magenta triangles) and right-handed selection 
mass (blue circles) versus the ratio of the GUT scale value of the wino soft mass to the bino soft mass, 
M 2 /Ml, for Model 2. 


with the neutralino, even in the presence of large stau mixing. At the same time it becomes 
difficult to accommodate an increasing within the constraints from {g — 2)^ so that 
no additional points with larger M 2 /Mi can be found in the 2a region. 

Conversely, for M 2 /Mi <C 1 the LSP becomes wino-like and drops down to very 
small values. As was explained at the end of Sec. 2, wino-like charginos and neutralinos are 
highly degenerate, which makes their detection very challenging at the LHC. Additionally 
light wino-like neutralinos require an extra dark matter component or production mech¬ 
anisms beyond those in the MSSM to satisfy For these reasons, we limit ourselves 
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here to the analysis of the parameter space yielding predominantly a bino-like LSP. 

Model 3 introduces a small difference between the right- and left-handed soft sfermion 
masses at the GUT-scale. The prototype we have in mind is a supergravity-based, GUT- 
scale model characterized by S'O(IO) boundary conditions [104, 105], where we assume a 
small positive H-term contribution, M|), from the extra (broken) U{1) [106, 107]. 

As usual, the GUT-defined soft squark masses (rug, rn^, rnjj), slepton masses {m\, 
m?^), and Higgs doublets masses can be parametrized in terms of the universal scalar 
mass in the fermionic 16 representation, mig, the universal scalar mass in the bosonic 10 
representation, ruio, and the H-term, M|,, so that 

2 _ 2 _ 2 _ 2 1 71 /r‘2 

~ = ”^16 + 

m% = = rriiQ — 3M^ 

^ ”^10 ^ 2M|, . (3.1) 

We scan mig, mio, and Mq in the ranges given in Table 2 and we neglect here the 
effects of including the right-handed sneutrino masses, as they have only slight impact on 
the low-energy spectrum [108]. Moreover, we scan tan/3 in the range 2 — 62 to increase 
the number of solutions, thus ignoring the requirement of successful 50(10) unification of 
the Yukawa couplings. A study of the impact of {g — 2)^ in 50(10) models including the 
Yukawa unification constraint can be found, e.g., in Ref. [109]. 

As was the case in Model 1, we assume GUT-scale unification of Mi and M2 and we 
leave M3 free to float, so to fit the Higgs mass without affecting the parameters that enter 
{g — 2)^, and so that at the EW scale a large, gluino-driven, sbottom mass can easily evade 
the bounds from direct searches at the LHG. 

Obviously, Model 1 is a subset of Model 3, so that the parameter space shown in 
Fig. 4(a) is common to both models. On the other hand, the extra freedom that comes in 
Model 3 from the right/left splitting at the GUT scale opens up an additional region of the 
low-scale parameter space, where the relic density is satisfied thanks to coannihilation of a 
bino-like neutralino and an almost degenerate right-handed slepton of the first or second 
generation. 

We show in Fig. 5(a) the physical and rrif^ distributions versus the neutralino 

mass for the right slepton-coannihilation region. As was the case in Model 1, an upper 
bound on can be derived, < 450 GeV, beyond which coannihilation becomes no 
longer efficient and starts to rise. The main contribution to 6 {g — 2)^ come from 

Eq. (2.5) and in this region /i ~ 4 — lOTeV. Figure 5(a) also shows that in this region the 
left-handed sleptons and, for many points, even the lightest staus, are not much heavier 
than the neutralino (and the lightest chargino, which is wino-like with ~ 2 m^). As we 
shall see in Sec. 4, this has important consequences when it comes to the LHG signatures. 

Model 4 generalizes Model 3 but extends the allowed GUT-scale ranges for the right- 
and left-handed sleptons. The prototype we have in mind is a supergravity model with 
minimal 5C/(5) boundary conditions at the GUT scale [110]. The GUT-defined soft squark 
and slepton masses are parametrized in terms of the common scalar soft mass for the fields 
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Figure 5: (a) The distribution of the physical left-handed selection mass (magenta triangles), right- 
handed selection mass (blue circles), and lightest stau mass (cyan diamonds) versus the neutralino mass for 
the points of Model 3 that are not in common with Model 1. All points satisfy the constraints of Table 3 
at the 2 (t level, (b) The distribution of the slepton physical masses in Model 4. 


belonging to the 5 representation, ms, and the common scalar mass in the 10, mig. The 
Higgs doublets’ soft masses are free, so that one has overall, 

mQ = 'rnu = rn% = ml^ 

ml) = ml = ml, (3.2) 

and ml ^^, m'jj^ are independent free parameters. Again, we assume that at the GUT scale 
Ml = M 2 = mi /2 and M3 is free-floating. The parameter ranges for ms, mio, m'j^^, m^j^ 
can be found in Table 2. 

Additionally, we also perform a scan of the parameter space of Model 4 after trading the 
GUT-scale inputs m'jj^ and m|^^ for the low-scale defined jjL parameter and the pseudoscalar 
pole mass, mA- We restrict the range of // to 10 — 2000 GeV, and at the same time we 
restrict the range of M3 to 500 — 2000 GeV (this scan is called Model 4-zoom in Table 2). 
We do this to counterbalance the tendency of our scans to find solutions characterized by 
very large values of fi and M3. In the spirit of phenomenology we will not consider issues 
of EW naturalness in this study. However, it is interesting to see if solutions more natural 
than the ones generated in Models 1-3 are possible. The scan of Model 4-zoom is designed 
to expose the region characterized by a Barbieri-Giudice measure [111, 112] roughly less 
than 1000 (see, e.g.. Ref. [93] for a discussion). 

Models 1 and 3 are sub-cases of Model 4, so that they share common regions of the 
parameter space. However, Model 4 introduces more freedom to the parameter space, with 
the consequence that there are regions for which the relic density is satisfied thanks to 
coannihilation of a left-handed slepton of the hrst or second generation and the neutralino. 
This is shown on the right of Fig. 5(b), where we plot the physical slepton mass distributions 
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for Models 4 and 4-zoom. Additionally, one can see on the left of Fig. 5(b), for ~ 
60 GeV, some solutions belonging to the /i-resonance region, which are efficiently explored 
by scanning the parameters in the limited ranges of Model 4-zoom. 

We want to point out here that the SU{5) boundary conditions (3.2) of Model 4 are 
not the ones exclusively giving the described low energy phenomenology. Because the 
GUT-scale value of M3 is allowed to be large and drive the squarks to decoupled values, 
the same solutions shown in Fig. 5(b) also apply to cases with different symmetries at the 
GUT scale, as would be the case, e.g., of Pati-Salam [113] boundary conditions. 

4 Constraints and projections for the LHC 

In this section we confront the GUT-defined scenarios discussed in Sec. 3 with the bounds 
from direct SUSY searches at the LHG. In Sec. 4.1 we briefly describe the numerical 
methodology we employ to derive the bounds from the 8 TeV run and we show our projected 
sensitivities for the 14 TeV run in SMS scenarios. In the following subsections we move on 
to discuss the impact of the implemented searches on the allowed parameter space of the 
considered models. 

4.1 Numerical implementation of LHC searches 

We numerically reproduce three LHG searches designed to explore the EW sector of the 
MSSM: the searches for EW chargino and neutralino production with 3 leptons in the 
final state by ATLAS and CMS [46, 47] (collectively called “3-lepton” hereafter), and a 
search for direct slepton pair production, sneutrino pair production, and slepton/sneutrino 
production with two opposite-sign leptons in the final state, by ATLAS [48] (dubbed as 
“2-lepton” hereafter). In 3-lepton analyses, CMS and ATLAS both reported at the end of 
the 8 TeV run small excesses in the observed events in different signal regions. Thus, each 
collaboration presented exclusion bounds slightly weaker than the expected ones, albeit not 
in the same signal regions. To take advantage of the stronger limits from each collaboration 
we adopt here a “best of’ strategy when we impose the limits from 3 lepton searches on 
the parameter space of our models. 

To recast the 3-lepton and the 2-lepton ATLAS searches, we employ the publicly avail¬ 
able code CheckMATE [114-122]. The analysis implemented in the package have been vali¬ 
dated by the code’s authors. We double checked by comparing the limits produced by the 
code with the official ones in three simplified models: chargino-neutralino production with 
VFZ-mediated decay into leptons; chargino-neutralino production with slepton-mediated 
decay into leptons; and left-handed slepton pair production. In all cases we found excellent 
agreement with the published results. 

The GMS 3-lepton search [47] is recast using the code designed by some of us and 
described in detail in Refs. [23, 36, 123]. For every point in the considered parameter space 
a set of 10® events is generated at the parton level with PYTHIAS [124], and the hadronization 
products are passed to the fast detector simulator DELPHES 3 [120] to reconstruct the 
physical objects. The GMS detector card is used, with the settings adjusted to those 
recommended by the experimental collaboration. Two kinematical variables proper of the 
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3-lepton search, invariant mass Mu and transverse mass Mt, are then constructed and 
used to divide the signal events into exclusive kinematical bins defined in the experimental 
paper. Finally, the acceptances/efficiencies are calculated as the fraction of all generated 
events that pass the applied cuts. The number of signal events is calculated as the product 
of the efficiency, luminosity and cross-section, where we use the NLO+NLL cross-sections 
provided by the LHC SUSY Cross Section Working Group [125]. 

The exclusion bounds are set according to a marginalized likelihood ratio method, 
with the Poisson likelihood distribution marginalized over the background and normalized 
to the background-only hypothesis as described in Sec. 3 of [123]. The total likelihood 
^LHC is given by the product of the likelihoods from each signal bin. 

Note, incidentally, that the ATLAS 2-lepton search uses instead signal regions that 
are not mutually exclusive. In that case, only the likelihood from the signal bin with the 
largest expected sensitivity is taken. As was mentioned above, when we combine the limits 
of different 3-lepton searches we assume that a point is excluded at the 95% C.L. when 
—21og£LHC >5.99 for at least one of the searches. 

We present in Fig. 6 the validation of the limits obtained with our code by comparing 
them to the official CMS 95% C.L. bounds. In Fig. 6(a) we show the SMS case of chargino- 
neutralino production with slepton-mediated decays in the “flavor-democratic” scenario, 
while in Fig. 6(b) we show the results for the case with no light sleptons. The points 
excluded by the likelihood function at the 99.7% C.L. are shown as gray dots, those excluded 
at the 95.0% C.L. as cyan circles, and those excluded at the 68.3% C.L. as blue triangles. 
Red squares indicate all other points. The solid black lines show the official CMS 95% C.L. 
exclusion limits from [47]. 

Besides the present LHC bounds we also calculate the sensitivity of the ATLAS 2- 
lepton and CMS 3-lepton searches at the LHC 14TeV run. We assume L = 300 fb~^ 
integrated luminosity. In each case we simulate the dominant SM backgrounds. For the 
3-leptons search these are WZ and tt production, as well as rare SM processes such as 
tiZ/W/H and triboson production. For the 2-lepton search the dominant backgrounds 
come from diboson production and tt production. Background events are generated at 
the LO using MadGraph5_aMC@NLD [126] and showered using PYTHIAS. The cross-sections 
are calculated at the NLO using MadGraph5_aMC@NL0. We generated 1.5 x 10® tt events, 
2 X 10® WZ events, and 3 x 10® rare SM process events for the 3-lepton search. For 
the 2-lepton search we additionally simulated 10® W~^W~ events and 10® ZZ events. 
The efficiencies for the background samples are derived applying the same experimental 
cuts used for the 8TeV run, and the number of background events is calculated as the 
product of the efficiency, luminosity and cross-section. The uncertainty in the number of 
background events is evaluated as the sum in quadrature of two terms: the uncertainty of 
the cross-section determination, given by MadGraph; and the statistical uncertainty of the 
efficiency determination with the Monte Carlo simulation.^ Note that this approach leads 
to a conservative estimate of future sensitivity, as one expects the eventual uncertainties 

®The uncertainty of the efficiency e is defined as ^ , where N is the total number of events 

generated in a Monte Carlo simulation. If the efficiency is equal zero, is reduced to = 1/N. 
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Figure 6: Exclusion confidence levels given by our 3-lepton likelihood function compared to the official 
CMS 95% C.L. exclusion bound [47] shown as a black solid line, (a) SMS case of chargino-neutralino 
production with slepton-mediated decays in the “flavor-democratic” scenario, (b) SMS case with with no 
intermediate sleptons decays. Gray dots are excluded at the 99.7% C.L., cyan circles at the 95.0% C.L., 
and blue triangles at the 68.3% C.L. Red squares indicate all other points. 


determined by the experimental collaborations using data-driven methods to be several 
times smaller. 

The background generation is first validated at 8 TeV by comparing the Monte Carlo 
results to the number of expected background events given in the experimental papers. In 
the course of the validation procedure the exclusion bounds for the “flavor democratic” and 
left-handed slepton SMS were also rederived using our background determination. In both 
cases we obtained good agreement of our procedure with the official experimental results. 

The projected exclusion bounds at 14 TeV are obtained by setting the number of 
observed events equal to the expected number of background events. In the ATLAS 2- 
lepton search discrimination between the signal and background yields is performed by 
means of the kinematical variable mT 2 [115, 121], with largest values of mT 2 probing large 
mass splittings between the slepton and neutralino. In the 8 TeV analysis the largest mT 2 
considered is mT 2 > 150 GeV. However, when the mass of the slepton increases, the 
mT 2 distribution for the signal falls more slowly than for the background well beyond this 
cut. Since at 14 TeV we expect to probe much larger slepton masses than in the 8 TeV 
case, we enhance the signal region by adding two new inclusive bins, mT 2 > 260 GeV and 
mT 2 > 310 GeV, to increases the sensitivity of the search in the high mass region. 

In Fig. 7(a) we present the projected 14 TeV sensitivity of the 3-leptons GMS search 
for chargino-neutralino production with slepton-mediated decays in the “flavor-democratic” 
scenario, with integrated luminosity L = 300 fb~^. The color code is the same as in Fig. 6. 
We calculate that the lower bound on the chargino mass for a neutralino LSP lighter than 
~ 900 GeV can be extended up to ~ 1400 GeV, which is a factor of two increase with respect 
to the 8 TeV result. In Fig. 7(b) we show the sensitivity of the ATLAS 2-lepton search for 
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Figure 7: Projected exclusion bounds for the LHC 14TeV run with 300 fb“^. (a) 3-lepton search for 
xfx 2 production with intermediate sleptons in the “flavor-democratic” SMS at CMS. (b) 2-lepton search 
for clS-l production at ATLAS. The color code is the same as in Fig. 6. The bounds from the 8TeV run 
are shown as dotted black lines. 


left-handed slepton pair production SMS. The difference in shape between the 8TeV and 
14TeV limits is due to the extra bins at large mT 2 , clearly indicating the importance of 
these bins in adding sensitivity in the heavy slepton region. 

4.2 Limits on GUT scenarios from the LHC 8 TeV run 

We show here the present LHC bounds on the parameter space of the models of Sec. 3. We 
obtain them by applying the simulation of the searches described in Sec. 4.1 to the model 
scans. For each point we simulate 10® events at LO using PYTHIAS for each of the relevant 
production mechanisms: X 2 X 1 and xf 2 X 3 4 for the 3-lepton searches and additionally 

and Xi'xr for the 2-lepton search. Cross sections are calculated at NLO 
using MadGraph5_aMC@NL0. 

We show in Fig. 8(a) the LHC bounds on the CMSSM-like scenario of Model 1. All 
points satisfy the constraints of Table 3 at the 95% C.L. The points excluded at the 
95% C.L. by the 3-lepton searches are shown as light gray triangles, those excluded by the 

2- leptons searches as gray circles, and those excluded simultaneously in both topologies as 
dark gray diamonds. The points that are presently not excluded by the LHC are shown 
as blue squares. The plot is presented in the (m, ±, mg^) plane, as the sensitivity of 

3- lepton searches to the chargino mass increases in this plane from right to left and is 
thus orthogonal to the sensitivity of 2-lepton searches for slepton pair production, which 
increases instead from top to bottom. Since Mi = M 2 = m,i /2 at the GUT scale, the 
approximate relation rriy^ ~ 0.5 m^± holds for all shown points. The exclusion criterion 
was described in Sec. 4.1. The surviving points have Ax^ = X^ — Xmin 5.99, where Xmin 
corresponds to the background only hypothesis. 


- 20 






























































































































BayesFITS (2015) 


BayesFITS (2015) 


Model 1 
=M2 

GUT unified sfermions 
LHC 8TeV 95% CL excl. 



m^± (GeV) 


O 500 

O 

I O 400 

£ 


Model 2 

Ml ^ M2 

GUT unified sfermions 
LHC 8TeV 95% CL excl. 




Excluded both 
Excluded 3 lep 
Excluded 2 lep 
Allowed 


m^± (GeV) 


(a) 


(b) 


BayesFITS (2015) 


> 

CD 

S 500 

E 


Model 3 
SO(10)-like 
LHC 8TeV 
95% CL excl. 





Excluded 3 lep 
Allowed 

rOg^ = m^± 
m, = 1/2 mY± 


600 700 

m + (GeV) 


800 



(c) 


(d) 


Figure 8: (a) The bounds from the LHC 8TeV run on the points of Model 1 in the (m ±, rrie,) plane. 
Points excluded at the 95% C.L. by one of the implemented 3-lepton searches are shown as light gray 
triangles, points excluded by the 2-lepton search are shown as gray circles, and points excluded in both 
topologies are shown as dark gray diamonds. Blue squares represent the points still unconstrained at the 
LHC. (b) Same as (a) for Model 2. (c) Same as (a) for the part of the parameter space of Model 3 that is 
not common to Model 1. (d) Same as (a) for Model 4. 


One can identify different areas of exclusion from left to right in Fig. 8(a). The few 
visible light gray triangles (and dark gray diamonds) on the left for m ± < 120 GeV are 
excluded by the 3-lepton search as xf and X 2 Sive rise to 3-body decays: e^xe^(or Ve) 
or ly^). We found that some models with > 120 GeV predominantly show 

large decay rates to on-shell staus, x^xf fiTfiiy-j-- As the chargino mass increases the 
search loses sensitivity due to the combined effect of the cross section drop and the fact 
that the mass splitting between x^ and fi is decreasing. 

The gray circles (and dark gray diamonds) for < 350 GeV are excluded by the 
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2-lepton bounds on direct left- and right-slepton production. The 2-lepton search loses 
sensitivity when the slepton is much heavier than the chargino and the dominant decay 
channel then becomes I —?■ yf r'z, which yields signihcantly softer final state leptons than 
those from I —)■ Ix- On the other hand, an alternative channel opens up in this region: 
Di —)• xf Since the production cross section of a sneutrino is comparable to that of the 
corresponding slepton, many points are within the sensitivity of the 2-lepton search. As a 
result all the parameter space with mg^ ~ < 250 GeV is excluded. 

Finally, the region on the right of the plot, for m ± > mgr, is excluded by the 3-lepton 

Xl 

searches, which are very efficient in detecting the X 2 xf U topology, unless the 

spectrum becomes excessively compressed in m ± and mg,. 

In Fig. 8(b) we show the equivalent exclusion plot for Model 2, which is characterized 
by M 2 7^ Ml at the GUT scale. The LHC bounds in the {m^±, mg^) plane do not show 
great differences with Model 1, with the exception of the region of the parameter space 
corresponding to m ± in the range 450 — 700 GeV, which is almost completely excluded 

Xl 

by the 3-lepton searches. As a matter of fact, at the origin of the more elongated shape of 
Model 2’s point distribution is the fact that the GUT-scale value of M 2 is allowed to be 
larger than Mi, giving rise to heavier charginos in the stau-coannihilation region, as heavy 
as ~ 1 TeV for some points. Solectron and chargino masses are, however, less compressed 
than in Model 1, so that 3-lepton searches are more sensitive to m^± < 700 GeV for 
Model 2. 

We showed in Sec. 3 that, when one considers GUT symmetry patterns that allow 
for disunifying the left and right-handed soft slepton masses, the allowed parameter space 
opens up and more ways to obtain the correct relic density become viable. Model 3 rep¬ 
resents scenarios for which symmetry considerations allow for small deviations from the 
universal case as happens, for example, in models of 50(10) SUSY GUTs with a positive 
O-term mass. In this case, besides the stau coannihilation and bulk regions, one hnds 
right-handed slepton/neutralino coannihilation. 

The LHC bounds on this additional part of the parameter space are shown in Fig. 8(c). 
The exclusion limits are given in this region by the 3-lepton searches, as the 2-lepton search 
for direct slepton production is still not sensitive to > 350 GeV. On the other hand, 
this scenario provides optimal conditions for the 3-lepton search, with sleptons of both 
“chiralities” such that m^ < m^^ < mj^ < m^±, so that the exclusion bounds are quite 
strong. One can see that points with m ± < 750 GeV are excluded, with the exception 

Xl 

of a few for which the bound is by ~ 100 GeV weaker. Those are points characterized by 
rufj < mg^, so that the enhanced values of BR(xi’^ —)■ fii^r) and BR(x2 —t fir) have the 
effect of weakening the search’s sensitivity. 

We want to point out, however, that the GUT-scale boundary conditions are such that 
a light stau in the spectrum is possible only for moderate values of tan /3, tan j3 < 25. In our 
scan, tan/3 is left a free parameter, but one ought to remember that a successful 50(10) 
unihcation pattern is likely to require large values of tan /3 to achieve Yukawa coupling 
unihcation, even if one allows for the possibility of substantial threshold corrections (see, 
e.g., Ref. [127]). Thus, the points not excluded with m^± < 750 GeV are less theoretically 
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motivated. 

Finally, we show in Fig. 8(d) the LHC bounds on Model 4, which is consistent with 
the SU{5) pattern of unification and others that allow for similar low-scale behavior. The 
parameter space consistent with {g — 2)^ opens up significantly in this case and broader 
regions can evade the present bounds. Besides the points where the mechanism for the relic 
density is right-handed slepton/neutralino coannihilation, as it was for Model 3, which are 
subject to the same bounds from the 3-lepton searches as the ones shown in Fig. 8(c), one 
can also see a stripe of points where the left-handed sleptons are almost degenerate with 
the neutralino, for chargino masses in the range 200 — 950 GeV. Here the search limit is 
weakened to < 600 GeV. This is due to the fact that, while the left-handed sleptons 
are compressed with neutralinos, in the same region the right-handed sleptons also happen 
to be very close in mass to the charginos, thus strongly reducing the efficiencies. 

In the regions above the diagonal, for m ± > 200 GeV, the strength of the 3-lepton 
searches is weakened by the fact that The dominant decay channels become 

X 2 xt hxiW^x^ for m ± ~ 450 — 600GeV and j^xf firfiVr for m ± ~ 200 — 

Xi Xi 

400 GeV. 

Finally, for the solutions comprising the /i-resonance region at m + ~ 100 — 120 GeV, g 
is limited to values in the 300 — 800 GeV range and stau-mixing is reduced so that for many 
points < m^± < The dominant decay channel, X 2 xf Vt-Vt is invisible to 
3-lepton searches. Pair production of sneutrinos of the first two generations do, however, 
decay to lepton -|- chargino and provide signatures within the reach of the 2-lepton search, 
but the sensitivity is not at the moment high enough to cover the whole region. One finds 
that it must be ~ > 240 GeV for the points with ~ 100 — 120 GeV. 

4.3 Projections for the LHC 14 TeV run 

In this section we investigate to what extent the second run at LHC will be able to explore 
the parameter space of the GUT-scale models defined in Sec. 3. We do so by applying the 
likelihood function described at the end of Sec. 4.1 to the model points in all our scans. We 
remind the reader that we assume 14 TeV center-of-mass energy and a target luminosity 
L = 300 fb"^ 

In Fig. 9(a) we show the projected 95% C.L. bounds for Model 1. The color code 
is the same as in Fig. 8. Remarkably, the parameter space can be excluded in its near 
entirety, predominantly by the 3-lepton search. It is interesting to note, however, that the 
regions of the parameter space for which the sensitivity of 3-lepton searches will remain 
insufficient, will be covered in complementarity by 2-lepton searches. This is the case, for 
example, of the points marked by gray circles at m ± ~ 100 — 180 GeV: their dominant 
decay chain is X 2 xf tittiV't, with stau and neutralino masses within ~ 20 GeV from 
one another. The resulting tau is so soft that detection will probably be a challenge even 
at 14 TeV. For the same points, however, the low-energy spectrum shows sleptons of the 
first two generations light enough to be easily detected in the 2-lepton search for direct 
pair production. 
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Figure 9: (a) The expected reach of the LHC 14TeV run on the points of Model 1 in the (m ±, mg,) 

X1 

plane. The color code is the same as in Fig. 8 . (b) Same as (a) for Model 2. (c) The projections for Model 2 
in the (M 2 /M 1 , M 3 /M 1 ) plane, (d) The projections for Model 2 in the (m ±, m^) plane. 

Xl 


The only points in the plot that survive our simulation, shown as blue squares, lie in 
the compressed spectra region ~ ~ ~ GeV. They also remain beyond the 

95% C.L. reach of the 2-lepton search, although our simulation places them within 68% 
C.L. reach. We remind the reader that our treatment of the 14TeV SM backgrounds is con¬ 
servative. When a more precise background determination is provided by the experimental 
collaborations, this region may be tested entirely with an even lower luminosity. 

If the gaugino mass universality condition is relaxed, as is the case of Model 2, new 
possibilities of evading the LHC exclusion bounds appear. In Fig. 9(b) we show the pro¬ 
jected bounds for Model 2. The points outside of the LHC reach, indicated by blue squares, 
are characterized by a GUT-scale ratio M 2 /M 1 < 1, as can be seen in Fig. 9(c) where we 
show the LHC projection in the plane of the ratios (M2/M1, M3/M1). These points feature 
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chargino and neutralino masses quite close to each other, as shown in Fig. 9(d) where we 
show the plane. Note that the distribution in the m^) plane (which we 

do not present here) show very similar behavior, namely the surviving points lie also on 
the compressed region for 2 -lepton searches, making this combination of parameters very 
challenging even for the 14TeV run. 

Recent studies that have looked into possible strategies for a more comprehensive 
coverage of SUSY spectra compressed in the EW sector at the LHC (and possibly the 
ILC) can be found, e.g., in Refs. [128-132]. The projected reach with 300fb“^ from, e.g., 
Ref. [129] is m ± ~ 250 GeV for a pure higgsino, which seems to fall short of probing these 
scenarios entirely. 

Incidentally, Fig. 9(c) also shows that it is difficult to generate points characterized at 
the same time by M 2 /Mi < 1 and very large values of M 3 at the GUT scale, as two-loop 
effects due to the gluino mass in the RGEs make the lightest stau tachyonic at the low 
scale. Hence the reduced density of points for M 3 /Mi > 8 . 

Moving on to Model 3, the part of the parameter space due to right-handed slepton 
and neutralino coannihilation is going to be entirely probed by 3-lepton searches with 
~ 100 — 110 fb~^ of integrated luminosity. Figure 10(a) shows that no point survives the 
cuts in this region. 

More interesting is the case of Model 4. We show in Fig. 10(b) the projected 95% C.L. 
exclusion bounds for Model 4 in the (m ±, mg^) plane. The only part of the parameter 
space that remains unconstrained corresponds to points with > m ±, characterized 
by large branching fractions for fiT or X 2 depending on whether the lightest 

stau is lighter or heavier than xf (and x^)- 

The surviving points situated at mg^ < 600 GeV should be in the future tested by the 
2-lepton search even with L = 300 fb~^ in the likely event that the uncertainties in the 
background determination will be eventually smaller than our estimate. The points with 
mg^ > 600 GeV, obviously outside of the 2 -lepton reach, are nonetheless characterized by 
large stau masses, and consequently larger branching fractions to the x^xf hxi lU^Xi 
channel. They remain beyond the reach of the 3-lepton search at 300 fb~^, but should 
eventually be tested with 3000 fb“^ [133]. 

To summarize our results, we show in Table 4 the parameters and branching fractions of 
3 benchmark points that survive our 14 TeV projections. Each point belongs to a different 
model. 

5 Summary and conclusions 

In this paper we have taken at face value the possibility that the ~ 3 cj anomaly in the 
measurement of 6 {g — 2)^ hnds its origins in the MSSM. If that is the case, the discrepancy 
should be to some extent confirmed by the New Muon g-2 experiment at Fermilab in the 
next few years. 

Given the present limits on SUSY masses from direct measurements at the LHG and the 
measurement of the Higgs mass at rrih — 125 GeV, a SUSY spectrum in agreement with all 
experimental constraints should feature sparticles charged under color significantly heavier 
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Model 3 that is not common to Model 1. (b) The LHC 14TeV projections for Model 4. In both panels the 
color code is the same as in Fig. 8. 


than the ones only charged under the EW gauge groups, whose masses must be in general 
of the order of a few 100s GeV to be in agreement with the measured 5 (g — 2)^. This 
raises the question as to what extent spectra characterized by this kind of hierarchy are 
consistent with the limits from the first run at the LHC and how deeply they can be further 
probed by the oncoming 14 TeV run. 

In generic parametrizations of the MSSM the issue has been examined in the literature 
and a brief reminder of the present LHC bounds on the parameter space subject to the 
additional requirement that the DM constraints are well satisfied by the lightest neutralino 
is presented in the first part of this paper. It is shown that large fractions of the parameter 
space can easily evade the present LHC limits and also lie beyond future sensitivity. 

On the other hand, realistic SUSY scenarios often present additional constraints on 
the parameter space due to the mechanism of SUSY breaking or additional symmetries. In 
gravity-mediated, GUT-defined scenarios it has been known for a while that models with 
gaugino non-universality can be at the same time in agreement with (g — 2 )^, direct limits 
from the LHC on colored sparticles, and the Higgs mass measurement, if the GUT-scale 
value of the gluino soft mass M 3 is substantially larger than Mi and M 2 . 

In this paper we exhaustively confronted the above scenarios with the exclusion bounds 
from direct SUSY searches at the LHC, particularly when additional constraints that come 
from the relic density and B physics observables are also taken into account. We considered 
four types of GUT-scale models characterized by non-universal boundary conditions. In 
the first we assumed GUT-scale universality of all the sfermions soft masses, as well as of 
the gaugino mass parameters Mi and M 2 . The gluino mass parameter, on the other hand, 
was allowed to be much heavier to boost the squark masses up to the multi- TeV regime, as 
required by the measured value of the Higgs boson mass, while leaving the EW part of the 


- 26 - 


















Benchmark 

BMl 

BM2 

BM3 

Model 

Model 1 

Model 2 

Model 4 

mo 

406.7 

280.0 

~ 

mio 

~ 

- 

269.8 

ms 

- 

- 

416.1 

tan/3 

22.9 

8.0 

18.4 

Ml 

747.6 

556.4 

759.4 

M2 

= Ml 

352.3 

= Ml 

Ms 

3427 

2860 

4227 

^0 

500 

-9.3 

1572 


= ^0 

= ml 

-(5.07 X 10^) 

2 

mir 

tiu 

= ml 

= ml 

-(5.01 X 10^) 

"'x? 

304 

219 

306 

m^± PS m^o 

577 

243 

579 

meL 

574 

288 

639 

rriej, 

498 

354 

310 

rrif^ 

329 

233 

772 

BF(xf) 

92% r Ur 
7% Ur T 

1% Ul 1 

100% f Ur 

100% Xl W 

BF(x^) 

92% f r 

7% Ur Ur 

1% Ul Ul 

100% f r 

93% Xl h 

6% [r 1 

1% X? ^ 

BF(eL) 

100% X? e 

54% xf 

28% xO e 
18% X? e 

67% Xl e 
22% x^ e 

11% xf 

BF(ejj) 

100% X? e 

100% Xl e 

67% Xl e 

BF(he) 

100% Xl 

51% xf e 
25% X2 

24% X? 

70% Xl 

20% xf e 
10% X2 

(3 lepton) 

1.8 

0.1 

0.29 
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Table 4: The model parameters for the benchmark points. All dimensionful quantities are given in GeV 
or GeV^. Also shown are the masses and branching fractions of the relevant particles in the EW sector as 
well as the values from the LHC searches at 14 TeV. 


spectrum relatively light to accommodate the measurement of {g — 2)^. In this scenario 
the proper value of DM relic density is obtained through neutralino LSP coannihilation 
with the lightest stau. 

This feature persists even when the gaugino mass universality condition is relaxed at 
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the GUT scale, Mi 7 ^ M 2 , which is the property of the second model we investigated. 
In the third model, slight relaxation of the universality condition on the sfermion masses, 
as happens within the framework of ^©(lO) GUT scenarios, introduces the possibility 
of efficient coannihilation with the right-handed sleptons. Finally, when the left- and 
right-handed slepton soft masses at the GUT-scale are treated as independent parameters, 
as is the case for example in our fourth, S'C/(5)-inspired model, also left-handed slepton 
coannihilation and Higgs resonance annihilation mechanisms become available. 

To investigate the impact of the LHC searches on these models we have simulated 
two kinds of searches to explore the EW part of the spectrum and therefore look for 
experimental signatures with varying number of leptons in the final state. The existing 2- 
and 3-lepton searches at ATLAS and CMS were recast using the publicly available code 
CheckMATE and a similar tool developed by some of us. 

On the other hand, the just started run H of the LHC, with a target center-of-mass 
energy of 14TeV, is expected to probe the EW sector of the MSSM much more efficiently. 
Thus, we derived predicted sensitivities of the 2-lepton and 3-lepton searches assuming a 
target luminosity L = 300 fb~^ and performing a detailed SM background simulation. 

Our analysis provides strong limits from the 8 TeV run on the parameter space of 
GUT-constrained scenarios consistent with [g — 2)^. However, a large number of model 
points are shown to evade the limits, thus leaving ample room for the explanation of the 
{g — 2)^ anomaly within GUT-scale SUSY models. 

In this regard, we have shown in this work that the parameter space surviving the 
bounds from the 8 TeV run falls within the sensitivity of the 14TeV run with 300 fb~^ 
projected luminosity virtually in its entirety. This opens up the interesting possibility 
that, if the (^f — 2 )^ anomaly is real and will be confirmed by future dedicated experiments, 
explanation within a large class of well motivated SUSY models will give sure signatures 
at the LHC or, alternatively, these models will have to be abandoned as an explanation for 
the {g — 2 )^ anomaly. 

A few small regions proved to be difficult to test even at the end of the 14TeV run, 
and they should be given special attention. 

• All of these models present large fractions of the parameter space where the neutralino 
is bino-like and almost degenerate with the lightest stau. This channel is notoriously 
difficult to test, because of the soft nature of the produced taus. Particular combinations 
of the input parameters have been shown to additionally conspire to push the chargino 
outside the reach of direct pair production and degenerate with the left selectron, so that 
both 3-lepton and 2-lepton searches cannot be used as a handle. 

• In models with Mi 7 ^ M 2 at the GUT scale not many combinations of Mi and 
M 2 allow one to obtain the correct relic density, and those that do feature in general 
M 2 /Ml > 1, so that they fall within the reach of 3-lepton searches. However, some models 
can have M 2 /Mi < 1, driving the EW spectrum to be highly compressed even for a 
predominantly bino-like neutralino. 

• Models with boundary conditions consistent with SU{5) or Pati-Salam are charac¬ 
terized by larger freedom in the parameter space, making them closer to what happens in 
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generic parametrizations of the MSSM with bino-like DM. The relic abundance constraint 
often requires at least one among the left- or right-handed selectron/smuon to be light and 
degenerate with the neutralino. However, the situation can arise where the other one and 
the lightest stau are several 100s GeV heavier than the slepton coannihilating. In this case 
2-lepton searches will be ineffective, as the model presents a compressed spectrum, and at 
the same time the sensitivity reach of 3-lepton searches will be curbed by the absence of 
intermediate, fairly light sleptons. X 2 decays predominantly into the Higgs bosons, 
requiring a much larger integrated luminosity. 

Increasing luminosity, reducing the background uncertainties and combining multiple 
searches statistically can be the first step to entirely probe the remaining model points. On 
the other hand, it is not excluded that SUSY hides exactly in “pockets” of the parameter 
space that are particularly challenging experimentally. Therefore, the effort should also be 
put in developing new search strategies that would have power to test such elusive spectra. 
In particular, we believe that finding efficient ways to increase the LHC sensitivity to 
compressed tau final states is of utmost importance as these taus seem to be the curtain 
behind which some SUSY scenarios might be hiding. 
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